Does Luttinger liquid behaviour survive in an 
atomic wire on a surface? 



L K Dash and A J Fisher 

Department of Physics and Astronomy, University College London, Gower Street, 
London WCIE 6BT 

E-mail: louise.dash@ucl.ac.uk, andrew.fisher@ucl.ac.uk 

Abstract. We form a highly simplified model of an atomic wire on a surface by 
the coupling of two one-dimensional chains, one with electron-electron interactions 
to represent the wire and and one with no electron-electron interactions to represent 
the surface. We use exact diagonalization techniques to calculate the eigenstates and 
response functions of our model, in order to determine both the nature of the coupling 
and to what extent the coupling affects the Luttinger liquid properties we would expect 
in a purely one-dimensional system. We find that while there are indeed Luttinger 
liquid indicators present, some residual Fermi liquid characteristics remain. 



Submitted to: J. Phys.: Condens. Matter 

PACS numbers: 71.10.Pm, 73.21.Hb, 79.60.Jv, 81.07.Vb 



Does Luttinger liquid behaviour survive in an atomic wire on a surface? 



2 



1. Introduction 

A purely one-dimensional metal is not expected to retain the Fermi liquid properties of 
its three-dimensional counterpart, which is characterized by a one-to-one relationship of 
its excitations to those of a non-interacting Fermi system and an electron quasiparticle 
lifetime that diverges to infinity near the Fermi surface. Instead, any interaction between 
electrons in a one- dimensional metal is expected to lead to a non-zero scattering phase- 
shift and to destroy completely the quasi-electron nature of the excitations [||]. The 
resulting "Luttinger liquid" state has separate charge and spin excitations, and quite 
different transport properties from those of a Fermi liquid. The Luttinger liquid state 
is an adiabatic continuation of the ground state of a simple one- dimensional model, the 
Luttinger model 0, in much the same way that a Fermi-liquid state can be adiabatically 
derived from the ground state of noninteracting fermions. 

How could we observe this fascinating state of matter? We cannot fabricate 
purely one- dimensional systems in the laboratory, but we can produce systems that 
are quasi-one- dimensional, to a greater or a lesser degree. In several of these, searches 
for Luttinger liquid behaviour have been undertaken. These include charge-transfer 
salts such as TMTSF having almost one-dimensional bandstructures , semiconductor 
quantum wires 0], edge states in the quantum Hall effect (which are chiral systems 
along which electrons can only propagate in one direction) |Q H, carbon nanotubes 
M, and atomic-scale quantum wires fabricated by self-assembly M or novel SPM-based 



surface lithographies [T^, [T^]. The important question for a theorist is whether such a 
quasi-one-dimensional system, interacting with a three-dimensional environment, would 
be expected to show Luttinger liquid behaviour or not, and if so, for what range of 
parameters. This question was given particular timeliness by the recent angle-resolved 
photoemission results of Segovia et al which have shown that photoemission from 
chains of Au atoms on Si(lll) exhibits features which may be signatures of the spin 
and charge collective modes indicative of Luttinger liquid behaviour. 

In this paper we begin to construct a theoretical framework within which we can 
answer these questions. We consider the simplest possible model of such a system — 
one in which the atomic wire is represented by a single-chain spinless Luttinger model, 
with electron-electron interactions, on a lattice. The environment is represented only 
by a chain with no interactions but which is otherwise identical. Electrons can hop 
between the two chains by hybridization of the electron states. Such a simple system 
does not capture the full complexities of electron scattering with the environment, nor 
can it represent the magnetic behaviour of a Luttinger liquid. However it provides 
a relatively simple starting point, and can already be used to calculate one-particle 
spectral functions such as those probed by photoemission experiments. 

In the following sections we explain the details of our model and our numerical 
methods in Section |^, and then describe our results in Section We find that some 
aspects of the model display Luttinger liquid-like behaviour while others do not; we 
then summarize these different properties in Section 
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2. Theoretical methods 

2.1. The Luttinger model 

We base our numerical calculations on Haldane's Luttinger model operator method 
for a chain of spinless idealized atoms with periodic boundary conditions. The Luttinger 
model is characterized by a linearized dispersion relation with slope fp and a complete 
separation of the populations of left- and right-moving particles (as in figure [^) plus 
the inclusion of an infinity of negative energy electrons p|, |^. The second of these 
characteristics makes the model exactly solvable by Haldane's bosonization method, as 
the boson commutation relations are then exact rather than approximate. The total 
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P= +1 
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Figure 1. Schematic representation of the Luttinger modeL p — +1, —1 indicate the 
populations of left- and right-moving fermions respectively. 

Luttinger model Hamiltonian can be split into three parts [Q: 

H = Ho + H2 + Ha. (1) 
Hq represents the non-interacting part of the Hamiltonian: 

p,k 

where q indicates momentum and p (= ±1) the branch index for the left- (—1) and right- 
moving (-1-1) fermions. The operator cL creates a fermion of type p with momentum k 
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and : . . . : implies normal ordering. The interaction parts of the Hamiltonian, H2 and 
^^4, are given respectively by 

^2 = -^X1^2(g)p+(g)p-(-g), (3) 

^4 = ^ ^ 94{q)Pp{q)Pp{-q)- (4) 

p,g 

The H2 term represents forward scattering between the left- and right-moving fermion 
branches, while the H4 term represents forward scattering within a momentum branch. 
Note however that [H2,Hq\ 7^ 0, so H2 can modify the ground state by excitation of 
particle-hole pairs, whereas H4 commutes with ^0 and so cannot modify the ground 
state. The density operators p are defined in terms of the fermion operators as 

Ppiq) = ■ 4,fc+/p-^ •= {^l,k+gCp,k - 5g,o(4,fcCp,fc)o) • (5) 

k k 

The complete Hamiltonian can be diagonalized via a Bogoliubov transformation 
to yield the bosonized form of the Hamiltonian 

^ = Z E ^o(^) ^ Pp(l)pp(-l) ■ +^ K (iV+ + N.f + vj (iV+ - N.f] (6) 

P, 97^0 

<j q 

where the transformed density operators p are related to the originals by a phase 0^: 
Ppiq) = Ppiq) cosh0g + p-piq) sinh0g. (8) 

There is a characteristic frequency ujq = \ivF+ g^iq) /2n)'^ — {g2iq) /27t)'^\^^'^ /\q\ associated 
with the transformed bosons, while the quantum numbers = + and J = 
— represent respectively the sum of the number of electrons on the positive and 
negative branches, and the difference between them (analogous to current). 
The three velocities in the Hamiltonian are related as follows: 

VnVj = Vq, 

vn = ^ = ^^oe-'^ (9) 

vj = voKp = foe+^'^, 
and are also related to the non-interacting Fermi velocity by 

, 5-4 + 6-2 
Vn =vf^ 7: , 

2^ (10) 

vj = Vf^ . 

zvr 

The requirement that if be a diagonalized version of H is ensured by the relationship 
between the Bogoliubov transformation phase 0g and the interaction functions gi: 

g20, ^ f T^VF + 94iq) - 92iq) \ 
\7TVf + g^iq) + g2iq) J 
= KM- (12) 
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The spinless Luttinger hquid parameter Kp is then obtained by taking the hmit of Kp{q) 
as q tends to zero. In all of the above the limit g ^ is implied where q is not explicitly 
included. 

In addition there are two further parameters related to Kp-. for a spinless Luttinger 
liquid they are given by 



a 



1 



1 
2 

2a 



- 2 



(13) 
(14) 



The Fermi liquid corresponds to Kp = 1 and a = and so departures from these values 
can be used to "measure" the extent of non-Fermi liquid behaviour, a is an exponent 
which governs the power-law dependence of all single-particle properties (for example 
the density of states, which varies as N{E) ^ \E — -Ef|"), as well as other properties of 
the system, including the d.c. resistivity, which varies as p{T) ^ T^^^ [|l], [l2| . 



2.2. Coupling two chains 

We are now able to couple two Luttinger chains (labelled by superscript A and B) 
together by allowing hopping between adjacent points on each chain with matrix element 
tr. 



^coupled ^ffA^jjB f^f{x)4j^{x) + ^^^(x)^'; 



(15) 



x,p 



Other schemes of inter-chain hopping are of course possible, and we expect that they 
will generate similar results. We neglect any "drag effects" of inter-chain interactions 
that would be generated by interaction terms analogous to H2 or H4 involving electrons 
on both chains. We have the choice of using either a set of basis states generated by the 
diagonalized boson operators 6j or the non-diagonalized operators aj. We choose the 
latter, as although this has the consequence that the ground states for a single interacting 
chain no longer consist of the relevant zero-boson basis states, it is considerably more 
convenient computationally. ■?/'''" are fermion creation operators given in the bosonized 
form by 

1 



71' 



[x) is a boson field operator given by 



[x] 



{^)Np + ^J2e{pq) 



2tx 



1/2 



e^-^^a, 



9' 



(16) 



(17) 



where the subscripts p refers to the branch index (±1) and c to the chain index, f/p ^ is a 
ladder operator whose form ensures the anticommutation properties of the final fermion 
operator '?/'^(x) despite the commutation properties of the constituent boson operators 
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Figure 2. Representation of the two chains. Electron-electron interactions are 
included on just one of the chains. 



In order for Lfp^c to produce anticommuting field operators on different chains, it is 
necessary to introduce a further phase factor into its definition, analogous to Haldane's 
phase factor ({p, Np, N_p) for ensuring anticommutation between the branches of a single 
chain. The total ladder operator component of equation (p!6D thus takes the form 

Up,c = C(P, Np, N_p)C'{c, N„ iV_e)|iVp,c + 1, iV_p,,), (18) 

where the subscript c = ±1 is a chain index. The anticommutation properties originate 
in the phase factors which can be written as 

C.=p,c=(-1)(^^^-). (19) 



2.3. Computational details 

The choice of a specific form for the interactions (72(0') and g^lq) is arbitrary, as the 
Luttinger liquid remains completely solvable for any interaction that fulfils certain 
conditions 0. We choose a Gaussian form for the interaction and set g2{q) = Qiio) = 
2TTV{q) with 

V{q) = Iexp{-2q^/r) (20) 

as this has the advantage that it maintains the same form in both real and momentum 
space. Other forms, such as a screened Coulomb interaction, would also be possible. The 
parameter / can be varied to control the "strength" of the interaction, and r controls 
the range of the interactions. 

Our model is completed by setting the interaction strength on chain B to zero 
but retaining interactions on chain A. With the non-interacting chain representing the 
surface and the interacting chain the wire, we therefore form our highly simplified model 
of an atomic wire interacting with a substrate. 

The size of the computational Hilbert space is restricted by allowing only one boson 
in each mode. While this successfully truncates the Hilbert space to a manageable size it 
does not affect the low energy properties in which we are interested, as we have checked 
that only higher energy excitations involve basis states with more than one boson in 
each mode. We also neglect the presence of the infinity of negative energy electrons. 
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i.e. c|A^eiectrons = 0) = 0. Howevei, we are still strongly limited in the size of the system 
we can handle: the Hilbert space scales exponentially with both length of chains and 
number of electrons. The largest system for which we have been able to calculate with 
arbitrary numbers of electrons is one with length L = 6 (i.e. 12 sites in total). In order 
to obtain both the eigenstates of the coupled system and its correlation functions we 



utilize the Lanczos method [O, [T^ 



2.4- Spectral functions 

We define the spectral function p{q, u) as 

p(g, u)p,, = --53 [GJ,(A;f + q,u + fi)], (21) 

71 

with q = kp — k, and the retarded Green's function defined as the double Fourier 
transform of 

G^,(a;,t) = -i0(t){(4,,(x,t)^t,(o,o))} (22) 

where the subscripts p and c again indicate branch and chain indices. It is possible 
to derive an analytic form for the spectral function of a single Luttinger liquid, either 



with or without spin, as in references ||I5|, [16], IT^]. Whereas for a Fermi liquid we would 



expect p(0, u) to be a delta function at the Fermi energy, the situation for a Luttinger 
liquid is quite different. Spectral weight is repelled from the Fermi surface owing to 
the virtual particle- hole excitations generated by the inter-branch scattering term g2, 
resulting in a broadened peak. As q is increased the Fermi liquid spectral function will 
merely broaden like reflecting the finite lifetime of electrons away from Ep, but for a 
Luttinger liquid there is zero spectral weight within a range ±fog of the Fermi energy. 
In addition, the negative frequency contribution is suppressed exponentially with q, and 
for a continuum Luttinger liquid the positive frequency contributions have a power law 
dependence. For uj > this is of the form 9{uj — Voq){u — v^q)'^'^ and for u; < it is of 
the form 6{—uj — Voq){—u! — Voqy, with 7 given by equation ([M]) |l|]. 

However, for our more complex system it is not easy to obtain an analytic form 
and so we resort to computational methods. We calculate the Green's function 

G'^{k,k',uj) = {N\ck\uj-H + ENY\l,\N) + {N\cl L + H-SNV^klN) (23) 



using Haydock's tridiagonal Lanczos-based procedure [0. In order to ensure 



convergence, we put u — > u + irj, where rj is an imaginary component of the energy 



roughly equal to the level spacing of the system |T^. |A^) is the A^-electron ground 
state, with energy ejy. 

3. Results 

In this section we present our numerical results, concentrating on L = 6 (our largest 
system). The results for L = 4 are very similar, and this encourages us to believe that, 
despite the small size of our system, finite-size effects are not dominating the results. 
All results shown here were calculated with vp = 1. 
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3.1. Ground and excited states 

The ground and neutral first excited states are calculated for different values of the 
interaction strength I as detailed above. As we have chosen the basis states for 
the system to be those before the Bogoliubov transformation, we find that the non- 
interacting ground state (i.e. / = 0) consists entirely of zero-boson basis states, whereas 
all the interacting systems (J > 0) have contributions from basis states that include 
bosons on one or both chains, as these basis states are not individually eigenstates of the 
interacting system. The contribution to the ground states from basis states containing 
either no bosons on either chain, bosons on the chain with explicit interactions only (i.e. 
chain A) and bosons on both chains are shown in figure ^ 
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Figure 3. Contribution of bosonic basis states to the ground states of the L = 6, 
N = 2, t± — 0.1 and L = 6, N = 2, t± = 0.5 systems as a function of interaction 
strength /. See text for a full explanation. 



The corresponding data for the boson contributions to a single isolated Luttinger 
liquid chain with otherwise identical parameters is also shown in figure ^ This matches 
almost exactly the contribution for the system's interacting chain (chain A) and thus 
indicates we are in a regime where the the two chains can be said to be weakly coupled 
for both coupling strengths = 0.1 and 0.5. 

The energies of the ground and first excited states are plotted in figure ^. We see 
that the ground state energy has a region of roughly linear dependence on / in the 
range < / < 0.5 and thereafter is rather less dependent on /. The first excited state, 
however, is linear in / for / < 2, at which point it becomes completely independent of /. 
The excitation energy AE for the non-interacting system is equal to 2t±, indicating that 
for this size of system the first excitation is one between the bonding and anti-bonding 
bands with zero net momentum change. As we increase I however, the bonding and 
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Figure 4. Energies for the ground state and neutral first excited state for the L — 6, 
N = 2, t± =0.1 system as a function of interaction strength /. 



antibonding bands distort to such an extent that eventually, at / ~ 2.0, the lowest 
excitation becomes not one between bands with zero net momentum, but an intra-band 
one with net momentum change 2ttv/L. For a Luttinger liquid, we would expect the 
velocity v in this expression to be the Luttinger velocity vq (as in equations @, etc.) 
rather than the Fermi velocity wp, leading to a continued /-dependence for AE. However, 
it is clear from figure ^ that this is not the case for this system: the excitation energy for 
/ > 2.0 remains constant and equal to It^vyIL = 7r/3. (For this system vp = 1.) This 
can be explained as follows — as J is increased the single-chain energy levels on chain A 
increase (since the interaction is repulsive), while those on the non- interacting chain B 
remain constant. It therefore becomes energetically preferable for the electrons to move 
to the non-interacting chain B rather than remain on the (interacting) chain A, until 
the point is reached at which the system is predominantly populated by electrons on 
chain B. This chain is still Fermi liquid-like, thus shifting the properties of the system 
as a whole back from partially Luttinger liquid-like to Fermi liquid-like. 

This is confirmed by an examination of the location of the electrons for both the 
ground and first excited states. Figure || shows a plot for L = 6, = 2,t± = 0.1, 
showing the chain to which the electrons "belong", again based on an analysis of the 
contribution of each of the basis states forming the eigenstate. As we have only 2 
electrons in this system, there are three possible distributions for the electrons: both on 
chain A ("A-like"), both on chain B ("B-like") or one on each chain ("equal"). 

First consider the ground states. We see that for the non-interacting ground state 
the electrons are distributed equally between the chains, and once interactions are 
switched on we evolve continuously towards a regime where basis states with electrons 
on chain A are strongly suppressed in favour of states with both electrons on chain B. 

The situation for the excited states is however quite different. We start at J = 
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Figure 5. Contributions to tlie ground and first excited states for the L = 6, N = 2, 
tj^ = 0.1 system from electrons on each chain as a function of interaction strength I. 
See text for a full explanation. 



with a state in which, hke the ground state, the electrons are equally distributed 
between the chains. This state is, however, 4-fold degenerate and immediately even weak 
interactions are introduced, this degeneracy is lifted and, in what is now the lowest state, 
the contribution from basis states with both electrons on the same chain are strongly 
suppressed. The first excited state for small values of / may therefore be thought of 
as "covalent" in nature, as opposed to the "ionic" -like states that contribute to the 
non-interacting state where both electrons are on the same chain. As interactions are 
increased further, the contribution from the ionic-like states begins to return. However, 
once we hit the point (at / = 2.0) where the nature of the excitations changes from 
zero net momentum change to Ak = 2ttv-p/ L, there is once more a discontinuity in the 
plot and the electron contributions from this point match those of the ground state, 
indicating that there is no net interchain electron hopping involved in these excitations. 

3.2. Spectral functions 

Results from the calculation of the spectral function p(g = 0, uj) are presented in figure 
^. There are three sources of the broadening of the peaks away from the delta function 
one would expect from a Fermi liquid. The first of these is due to the finite value of the 
imaginary energy rj we have used in our calculations, leading to broadening of the peaks 
into Lorentzians with width ~ rj. The second results from the splitting into bonding and 
anti-bonding bands due to the inter-chain coupling and is equal to 2t±. The remaining 
broadening, however, is due to the removal of spectral weight from the Fermi energy at 
cij = 0, and is indicative of Luttinger hquid behaviour in the coupled system. 

We can also calculate p{q ^ 0, u) although for our extremely small system, we 
are severely restricted in available values of q. As we would expect for a Luttinger 
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Figure 6. Hole (left peak) and electron (right peak) contributions to the spectral 
function p{q ~ 0,0;). t± = 0.1, rj = 2t±. In order to show the splitting clearly, the 
electron and hole contributions arc shown as equal in magnitude. 



liquid, there is a strong suppression of the hole contribution for g > 0. Owing to the 
discrete nature of this system, it is not possible to calculate the Luttinger parameter Kp 
direct from the slopes of the electron and hole contributions as in references |T5|, IT6|, [17 



However, we can still derive a value for Kp from the separation of the electron and 
hole peaks as follows: For a non-interacting coupled chain system, the electron and hole 
contributions to p{q = 0,uj) are separated by 2t±, whereas for a single Luttinger liquid 
the contributions to p(g ^ 0, u) are separated by 2t>og. We define an effective Vq in 
terms of the peak splitting Au by the relation 

Auj = 2{vfq + t^). (24) 

Then, since we have g2{q) = g^iq), we can use Vq^ as a measure of vq and equate 
equations (H) and (|l^) to give an expression for the Luttinger parameter: 



k: 



-20 



-^0 ) 



,eff 



,efr2 



f F + 



Vf 



(25) 
(26) 



since fp = 1. The results of this calculation and the corresponding values of the 
parameter a (equation ([13|) ) are plotted in figure 0. 



4. Conclusions 



Our calculations approach the question of whether our model exhibits Luttinger liquid 
behaviour from several complementary directions. We have presented results for the 
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Figure 7. Luttinger parameters Kp and a for the coupled chain system, both L = 6 
and L = 4, and for an isolated continuum limit (L — > oo) chain. 



nature of the ground state, for the nature of excited states (charged and neutral) and 
the related response functions governing the response of the system to excitation. We 
have also estimated the exponents in certain power-laws characterizing the Luttinger 
liquid state. 

Our results display a certain ambiguity with regard to whether or not this system 
is a Luttinger liquid. Roughly speaking, the evidence from the bosonic nature of the 
ground state, the charged excited states, and the values of the scaling parameters is 
all broadly consistent with Luttinger liquid behaviour, modified by the presence of the 
second chain. 

Specifically, the one-particle spectral function is consistent with Luttinger liquid 
behaviour. At g = 0, the total spectral function is broadened by separation of the 
electron and hole parts — this separation is much greater than the 2t± expected in a non- 
interacting two-chain system. For q ^ 0, there is additional separation of the electron 
and hole parts, consistent with the further suppression of spectral weight near the Fermi 
energy. (For g > 0, corresponding to probing the occupancy of states with k > kp — 
see equation (^) — the hole part of the spectral function is strongly suppressed, even 
for relatively strong interactions.) Further support for the Luttinger liquid hypothesis 
comes from the composition of that portion of the ground state in which bosons are 
excited on chain A: this fraction closely tracks the fraction that would be expected in 
an isolated single-chain Luttinger liquid (Figure^). 

Moreover, we are able to derive indirectly from p{q > 0,uj) values for two Luttinger 
liquid parameters: Kp (which determines the ratio of the different velocities) and a 
(which determines the asymptotic behaviour of the correlation functions). Although 
our analysis presupposes the validity of relations characterizing the Luttinger liquid 
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state (equation (|25|)), the results are self-consistent in the sense that they show Kp < 1 
and a > 0, i.e., departures from the Fermi- liquid values. The degree of departure from 
Fermi-liquid values is, however, somewhat less than would be expected for a single-chain 
system with the same interaction strength (see Figure |^). 

However, we also have evidence that there is some residual Fermi liquid behaviour. 
This comes principally from the distribution of charge across the two chains, as a 
function of interaction strength, in the ground and neutral excited states. Transfer 
of charge to the non-interacting chain at higher values of I results in the resurrection 
of a Fermi-liquid like neutral excitation spectrum (although we are still able to extract 
non-Fermi liquid values for the Luttinger parameters for these interaction strengths, as 
explained above). 

This implies that while the neutral excitation spectrum undergoes a transformation 
to Fermi liquid behaviour at large J, the charged excitation spectrum does not. It 
remains to be seen how the apparent conflict between these different views on the same 
underlying quantum state is resolved. 

Finally, it should be noted that in our system the surface is represented as metallic, 
although in order to achieve true one-dimensional conduction in an experiment it is 
much more likely that a semiconducting or insulating substrate would be employed, as 
in the work of [||. One might expect that with a gapped substrate, the elimination 
of hybridization with substrate states at the Fermi energy would promote Luttinger 
liquid formation by comparison with the metallic case. So our finding that some aspects 
of Luttinger liquid spectral behaviour persist even in the metallic case may also have 
implications for recent and future experiments. We plan further theoretical work on 
semiconducting substrates. It would also be interesting to study the influence of wire- 
substrate interactions other than simple hybridization on Luttinger liquid formation. 
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